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We analyze the effect of primordial density perturbations on the cosmic QCD 
phase transition. According to our results hadron bubbles nucleate at the cold 
perturbations. We call this mechanism inhomogeneous nucleation. We find the 
typical distance between bubble centers to be a few meters. This exceeds the 
estimates from homogeneous nucleation by two orders of magnitude. The resulting 
baryon inhomogeneities may affect primordial nucleosynthesis. 
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The order of the QCD transition and the values of its parameters are 
still under debate. Nevertheless there are indications from lattice QCD cal- 
culations. Quenched QCD (no dynamical quarks) shows a first-order phase 
transition with a small latent heat, compared to the bag model, and a small 
surface tension, compared to dimensional arguments El. We assume that the 
QCD transition is of first order and that the values from quenched lattice QCD 
(scaled appropriately by the number of degrees of freedom) are typical for the 
physical QCD transition. Based on these values and homogeneous bubble nu- 
cleation a small supercooling, Age = 1 — Tf/IL ~ 10"'*, and a tiny bubble 
nucleation distance, d^uc ~ 1 cm, would follow □. The actual nucleation tem- 
perature is denoted by Tf , and the thermodynamic transition temperature by 
Tc « 150 MeV. 

We argue i that the assumption of homogeneous nucleation is violated 
in the early Universe by the inevitable density perturbations from inflation 
or from other seeds for structure formation. Those fluctuations in density 
and temperature have been measured by COBEa to have an amplitude 
6T/T 10~^. The effect of the QCD transition on density perturbations 
and gravitational wavesQ has been studied previously, while we investigate the 
effect of the density perturbations on the QCD phase transition here. 

First-order phase transitions normally proceed via nucleation of bubbles 
of the new phase. When the temperature is spatially uniform and no signif- 
icant impurities are present, the mechanism is homogeneous nucleation. The 
probability to nucleate a bubble of the new phase per time and volume is ap- 
proximated by r « exp[—S{T)]. The nucleation action S is the free energy 
difference of the system with and without the nucleating bubble, divided by 
the temperature. 

Nucleation is a very rapid process, compared with the extremely slow 
cooliiigjDf the Universe. The duration of the nucleation period, Atnuc, is found 
to belil 

^1/3 

(1) 



The time i/ is defined as the moment when the fraction of space where nu- 
cleations still continue equals 1/e. The heat flow preceding the deflagration 
fronts reheats the rest of the Universe. We denote by Whoat the effective speed 
by which released latent heat propagates in sufficient amounts to shut down 
nucleations. In practice, Vdef < ^^heat < Cg, pjihere Vdd is the velocity of the 
deflagration front and Cg is the sound speed E3. In the unlikely case of deto- 
nations Wheat should be replaced by the velocity of the phase boundary in all 
expressions that follow. 
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The mean distance between nucleation centers, measured immediately af- 
ter the transition completed, is 



^nuc,hom ^^hcat^^nuc- (2) 

This nucleation distance sets the spatial scale for baryon number inhonio- 
geneities. 

Lattice simulations imply that in real-world QCD the energy density must 

change very rapidly in a narrow temperature interval. This can be seen from 

1 11 

the microscopic sound speed in tlmquark phase, Cs = {dp/de)^ . Lattice QCD 
indicates that 3Cs(Tc) = 0(0.1)11^1. Thus, the cosmological time-temperature 
relation is strongly modified already before the nucleations, due to 

where tn = = (3A/p[/87r£q)^''^ with eq being the energy density in the 
quark phase. This behavior of the sound speed increases the nucleation dis- 
tance because of the proportionahty Ainuc (x l/[3Cs(T/)]0. 

In the thin-wall approximation the nucleation action has the following 
explicit expression: 

for small supercooling. Assuming further that Cg does not change very much 
during supercooling, the following relation holds for the supercooling and nu- 
cleation scales: 

A^sc ^sc ^ Q {K\ 

At " A ~ ^ ^ ' 

^^nuc ^nuc 

Here we denote by A a relative (dimensionless) temperature interval and by 
At a dimensionful time interval. S = S{Tf) is the critical nucleation action, 
5 = 0(100). 

Surface tension and latent heat are provided by lattice sirnulations with 
quenched QCD only, giving the values a = O.OIST^, I = lAT^K Scaling the 
latent heat for the physical QCD leads us to take I — 3T^. 

With these values for the latent heat and surface tension, the amount of 
supercooling is Age — 2.3x10"^. FromEq. (||) it follows that Anuc = 1.5x10"^. 
Substituting 3c^ = 0.1 into Eq. (|), we find Atnuc = 1-5 x lO'^t^ for the 
duration of the nucleation period. The nucleation distance depends on the 
unknown velocity Wheat in Eq. (||). With the value 0.1 for fhcat, the nucleation 
distance dnuc hom would have the value 2.9 x 10~^dH- One should take these 
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values with caution, due to large uncertainties in / and cr. As our reference set 
of parameters, we take: Age = 10^'', Anuc = 10^^, Ainuc = 10~^i//. 

In the real Universe the local temperature of the radiation fluid fluctu- 
ates. We decompose the local temperature T{t,x.) into the mean temperature 
T{t) and the perturbation (5T(t,x). The temperature contrast is denoted by 
A = 5T/T. On subhorizon scales in the radiation dominated epoch, each 
Fourier coefficient A(t, k) oscillates with constant amplitude, which we denote 
by Axik). Inflation predicts a Gaussian distribution, 

1 / 1 A2 \ 
p(A)dA = ^= exp — — ^ dA . (6) 

We find El for the COBE normalized i rms temperature fluctuation of the 
radiation fluid (not of cold dark matter) A™ = 1.0 x 10"'' for a primordial 
Harrison-Zel'dovich spectrum. The change of the equation of state prior to 
the QCD transition modifies the temperature-energy density relation, A = 
cl6e/{e + p). We may neglect the pressure p near the critical temperature 
since p £q at Tc. On the other hand the drop of the sound speed enhances 

the amplitude of the density fluctuations proportional to Cs 13. Putting all 
those effects together and allowing for a tilt in the power spectrum, the COBE 
normalized rms temperature fluctuation reads 

A^^- « 10-^3c^)3/4 (^^ j , (7) 

where ko — {aH)o- For a Harrison-Zel'dovich spectrum {n = 1) and Bc^ = 0.1, 
we find A™*^ w 2 x 10"^ 

A small scale cut-off in the spectrum of primordial temperature fluctua- 
tions comes from collisional damping by neutrinos lij'Q. The interaction rate 
of neutrinos is GpT^. This has to be compared with the angular frequency 
Cskph of the acoustic oscillations. At the QCD transition neutrinos travel freely 
on scales Zj, « 4 x 10~^dH- Fluctuations below the diffusion scale of neutrinos 
are washed out, 

rta 



diff 



Z^(t)dt 



TxlQ-'^dH. (8) 



In Ref.S the damping scale from collisional damping by neutrinos has been cal- 
culated to be kP^ = lO^H at T = 150 MeV. The estimate (|) is consistent with 
this damping scale. We assume /smooth — lO^^'dn. The compression timescale 
for a homogeneous volume ^smooth ~ ""'smooth/cs ~ 10~^tH- Since 

6t ^ Aiimc the temperature fluctuations are frozen with respect to the time 
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Figure 1: Sketch of a first-order QCD transition in the inhomogeneous Universe. At ti 
the first hadronic bubbles (H) nucleate at the coldest spots (light gray), while most of the 
Universe remains in the quark phase (Q). At t2 the bubbles inside the cold spots have merged 
and have grown to bubbles as large as the temperature fluctuation scale. At tg the transition 
is almost finished. The last quark droplets are found in the hottest spots (dark gray). 



scale of nucleations. As long as /smooth exceeds the Fermi scale homogeneous 
bubble nucleation applies within these small homogeneous vplumes. This is a 
crucial difference to the scenario of heterogeneous nucleation Eil, where bubbles 
nucleate at ad hoc impurities. 

Let us now investigate bubble nucleation in a Universe with spatially inho- 
mogeneous temperature distribution. Bubble nucleation effectively takes place 
while the temperature drops by the tiny amount Anuc- To determine the mech- 
anism of nucleation, we compare Anuc with the rms temperature fluctuation 

1. If A™ < Anuc, the probability to nucleate a bubble at a given time is 
homogeneous in space. This is the case of homogeneous nucleation. 

2. If A™ > Anuc, the probability to nucleate a bubble at a given time is 
inhomogeneous in space. We call this inhomogeneous nucleation. 

The quenched lattice QCD data and a COBE normalized flat spectrum 
lead to the values Anuc ~ 10"® and A™ ~ 10"^. We conclude that the 
cosmological QCD transition may proceed via inhomogeneous nucleation. A 
sketch of inhomogeneous nucleation is shown in Fig. y. The basic idea is 
that temperature inhomogeneities determine the location of bubble nucleation. 
Bubbles nucleate first in the cold regions. 

The temperature change at a given point is governed by the Hubble expan- 
sion and by the temperature fluctuations. For the fastest changing fluctuations, 
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with angular frequency //smooth, we find 




3c? + o(a.| 



) 



(9) 



The Hubble expansion is the dominant contribution, as typical values are 3c^ — 
0.1 from quenched lattice QCD and A™''iH/'^i ~ 0.01 from the discussion 
above. This means that the local temperature does never increase, except by 
the released latent heat during bubble growth. 

To gain some insight in the physics of inhomogeneous nucleation, let us first 
inspect a simplified case. We have some randomly distributed cold spheres of 
diameter ^smooth with equal and uniform temperature, which is by the amount 
A™Tc smaller than the again uniform temperature in the rest of the Universe. 
When the temperature in the cold spots has dropped to Tf, homogeneous 
nucleation takes place in them. Due to the Hubble expansion the rest of the 
Universe would need the time Atcooi = i_f/A™'^/3c^ to cool down to Tf. Inside 
each cold spot there is a large number of tiny hadron bubbles, assumed to grow 
as deflagrations. They merge within Atcooi if Anuc < (wdof/whGat)A™^. This 
condition should be clearly fulfilled for our reference set of parameters. Thus 
the cold spots have fully been transformed into the hadron phase while the 
rest of the Universe still is in the quark phase. The latent heat released in a 
cold spot propagates in all directions, which provides the length scale 



If the typical distance from the boundary of a cold spot to the boundary of a 
neighboring cold spot is less than Zhoat, then no hadronic bubbles can nucleate 
in the intervening space. In this case the nucleation process is totally domi- 
nated by the cold spots, and the average distance between their centers gives 
the spatial scale for the resulting inhomogeneities. In the following analysis 
for a more realistic scenario we concentrate in this case, Zheat > Smooth- 

The real Universe consists of smooth patches of typical linear size /smooth, 
their temperatures given by the distribution (^. As discussed above, the merg- 
ing of tiny bubbles within a cold spot can here be treated as an instantaneous 
process. The fraction of space that is not reheated by the released latent heat 
(and not transformed to hadron phase), is given at time t by 



where we neglect overlap and merging of heat fronts. At time t heat, coming 
from a cold spot which was transformed into hadron phase at time t' , occupies 



^heat — ^l^hcat A^cool- 



(10) 




(11) 
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the volume V{t,t') — (47r/3)[/smooth/2 + Whcat(t — The other factor in 

Eq. (pi]), Fihri, is the volume fraction converted into the new phase, per physical 
time and volume as a function of the mean temperature T = T{t). Fihn is 
proportional to the fraction of space for which temperature is in the interval 
[T[, Tf{l + dA)] . This fraction of space is given by Eq. (|) with A = Tf /T - 1. 
Rewriting dA by means of Eq. (^) leads to the expression 

Tihn = 3c^^— ^ piA = ^ - 1), (12) 

J IH V'smooth 

where the relevant physical volume is Vsmooth = (47r/3)(Zsmooth/2)'^. 

The end of the nucleation period, ^ihn, is defined through the condition 
/(tihn) = 0. We introduce the variables iV = (1 - Tf/T)/A'r^'' and Af = 
N{titia)- Since Cs may be assumed to be constant during the tiny temperature 
interval where nucleations actually take place, we find from Eq. (^): 1 — t/tihn ~ 
2 / (Sc^) A™^ {N — J\f). Putting everything together we determine Af from 

Jhea^ / aN^^l^-^p^+N-M) =1. (13) 
'smooth V 27r V 'heat / 

The expression /hcat/^smooth = 2whoat(3c^)"^/''(fc/fco)'""^^/^ is valid for the 
COBE normalized spectrum. For ^hcat/4mooth = 1,2,5,10 we find Af ~ 
0.8, 1.4, 2.1, 2.6, respectively. 

The effective nucleation distance in inhomogeneous nucleation is defined 
from the number density of those cold spots that acted as nucleation centers, 
rfnucihn = We find 



ihn 



tihn 



Fihn(i)di 







-1/3 



(14) 



= [3 (1 _ erf(AA/V2)]-i/3 l,^,,,^,. (15) 

TT 

With the above values ^hcat /^smooth = 1,2,5, 10 we get dnuc,iim = 1-4, 1.8,3.0, 
4.8 X ^smooth, where /smooth « 1 m. 

For a COBE normalized spectrum without any tilt and with a tilt of n — 1 = 
0.2 (where (fcsmooth/fco)*''^ ~ 25), together with 3cl — 0.1 and Whoat = 0.1, we 
find the estimate /hcat/'smooth ~ 0.4 and 9, correspondingly. Notice that the 
values of Wheat and 3Cj are in principle unknown. Anyway, we can conclude 
that the case /heat > ^smooth is a realistic possibility. 

With 2uhcat(3c^)~^/'*(10~'*(iH//smooth) < 1 and without positive tilt we are 
in the region Ihcat < ^smooth, where the geometry is more complicated and the 
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above quantitative analysis does not apply. In this situation nuclcations take 
place in the most common cold spots {Af ~ 1), which are very close to each 
other. We expect a structure of interconnected baryon-depleted and baryon- 
enriched layers with typical surface Zgj^ooth ^^'^ thickness Idcf = Wdcf Aicooi- In 
between dnuc,hom would be the relevant length scale of inhomogeneities. An 
accurate analysis of this case requires computer simulations, which is beyond 
the scope of the present work. However, it is clear that the result will be 
different compared with homogeneous nucleation. 

We emphasize that inhomogeneous and heterogeneous nucleation are 
genuinely different mechanisms, although they give the same typical scale of a 
few meters by chance. If latent heat and surface tension of QCD would turn 
out to reduce Age to, e.g., 10~^, instead of lO""*, the maximal heterogeneous 
nucleation distance would fall to the centimeter scale, whereas on the distance 
in inhomogeneous nucleation this would have no effect. 

We have shown that inhomogeneous nucleation during the QCD transi- 
tion can give rise to an inhomogeneity scale exceeding the proton diffusion 
scale (2 m at 150 MeV). The resulting baryon inhomogeneities could provide 
inhomogeneous initial conditions for nucleosynthesis. Observable deviations 
from the element abundances predicted by homogeneous nucleosynthesis seem 
to be possible in that caseEj. 

In conclusion, we found that inhomogeneous nucleation leads to nucleation 
distances that exceed by two orders of magnitude estimates based on homoge- 
neous nucleation. We emphasize that this new effect appears for the (today) 
most probable range of cosmological and QCD parameters. 
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